Spin-Based Quantum Computers made by Chemistry: Hows and Whys 
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This introductory review discusses tlie main problems facing tiie attempt to build quantum in- 
formation processing systems (like quantum computers) from spin-based qubits. We emphasize 
'bottom-up' attempts using methods from chemistry. The essentials of quantum computing are 
explained, along with a description of the qubits and their interactions in terms of physical spin 
qubits. The main problem to be overcome in this whole field is decoherence - it must be considered 
in any design for qubits. We give an overview of how decoherence works, and then describe some 
of the practical ways to suppress contributions to decoherence from spin bath and oscillator bath 
environments, and from dipolar interactions. Dipolar interactions create special problems of their 
own because of their long range. Finally, taking into account the problems raised by decoherence, 
by dipolar interactions, and by architectural constraints, we discuss various strategies for making 
chemistry-based spin qubits, using both magnetic molecules and magnetic ions. 



I. INTRODUCTION 

The basic question addressed in this short review is: 
how might one build a spin-based quantum computer us- 
ing the 'bottom-up' methods of chemistry? Before even 
starting on this topic, the reader might suppose that we 
should first answer two other questions, viz. (i) can one 
build a quantum computer at all? (ii) why would this be 
a subject for chemists? 

The answer to the first question is not known. Al- 
though one can easily imagine various schemes for do- 
ing quantum information processing, to actually build 
a 'Quantum Information Processing System" (QIPS) 
one has to deal with a really quite fundamental prob- 
lem in Nature, usually called the 'decoherence problem'. 
Quantum computing involves manipulating 'multiply en- 
tangled wave- functions', involving many different sub- 
systems in the QIPS - in these terribly complex quan- 
tum states, distributed over the whole QIPS, the indi- 
vidual wave-functions of the individual components lose 
all meaning and only the wave function of the entire sys- 
tem is physically meaningful. The problem is that such 
states are extremely delicate, and can be destroyed by 
very weak interactions with their surroundings - and yet 
in order to use them we need to be able to probe and 
manipulate them. Physicists and chemists working in 
'nanoscience' have come up with many designs, and there 
has been considerable success in building QIPS involving 
just a few entangled 'qubits' (the simplest kind of sub- 
system, involving only two quantum states). This success 
has simply made the importance of the decoherence prob- 
lem more obvious - solving it will be the key that unlocks 
the QIPS door (and doubtless many other doors as well). 

Concerning the second question - so far the most 
successful designs for QIPS have been made from ion 
traps 1], but it is not yet clear how one might scale 
these up to the larger systems, involving many entan- 
gled qubits, that are needed for quantum computation. 
The most likely designs may be hybrids between 'solid- 
state' circuits of qubits, with communication effected by 



photons - nobody really knows yet. There are then 2 
good reasons for concentrating on spin-based systems, 
made using bottom-up methods. The first is that a de- 
vice which only uses spin degrees of freedom avoids mov- 
ing charges around - this is good because moving charges 
invariably causes decoherence. The second is that if 
the basic sub-units are made by Nature (in the form of 
molecules or ions) then quantum mechanics guarantees 
they will all be identical, provided we can eliminate im- 
purities, defects, etc. This of course takes us into the 
realm of chemistry. 

In writing this review we have assumed a reader who 
may not be too familiar with this subject, but who is in- 
terested in understanding the main ideas and problems, 
and who is also looking for suggestions for future work in 
this area. We therefore begin, in section 2, by describing 
briefly what a QIPS is, and how one would describe its 
physical components if they were indeed made from spin 
qubits using chemistry. This description also includes 
a discussion of the physical mechanisms responsible for 
decoherence in such systems. Given the importance of 
decoherence, we have devoted all of section 3 to how it 
works. The emphasis is on simple pictures, which often 
convey a better feeling for what is going on than com- 
plex mathematics. Then, in section 4, we outline some 
features of the chemistry, and a few desiderata for any 
future architecture. This allows us to make some sug- 
gestions for future work. This last section is necessarily 
open-ended - at present nobody knows which designs may 
ultimately succeed. 



II. DESCRIBING THE SYSTEM 

Physicists and chemists on the one hand, and com- 
puter scientists on the other, have quite different ways of 
describing a system whose ultimate aim is to do comput- 
ing (quantum or classical). Each kind of description is 
valuable, but the differences often cause misunderstand- 
ing. In this section we describe and then connect the two 
kinds of description; we then discuss an effective Hamil- 
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FIG. 1: Cartoon view of a Turing machine: a head (shown 
shaded) on a control loop reads and writes onto a linear 
"tape". The states of both are binary (shown here as spins). 
We do not show the register or the program parts of the ma- 
chine. 



tonian for a set of qubits made from real materials. 



A. Description in Computational Language 

Computer scientists like to reduce physical computa- 
tional systems to " Turing machines" [3| ; a clear physical 
discussion of this was given by Feynman[3|, |j]. In a Tur- 
ing machine, a "tape" moves forward or backwards past a 
"head" ; the head has the job of writing in or reading out 
particular discrete states (labelled by different symbols) 
on the tape (and it may also erase states) . There is also a 
table of instructions which tell the tape and head what to 
do (the "program"), and a "stable register" or memory, 
which stores the state of the machine at any given time. 
(Turing thought of this as being analogous to the state 
of mind of a sentient being). The net result is shown in 
Fig. [1] Note that the state of a Turing machine at any 
time can be understood in binary code as a string of I's 
and O's (or a sequence of on/off states, or up/down, or 
heads/tails, etc.). 

Even though real computers may not be constructed 
like a Turing machine (indeed Turing machines provide 
a very clumsy model for a realistic computing archi- 
tecture), nevertheless any computer can be mapped, in 
its logical and operational structure, to a Turing ma- 
chine. Moreover, Turing showed that a rather simply 
constructed machine of this kind (since made even sim- 
pler) could do any discrete computation. Thus all digital 
computers can be described in terms of such "universal 
(classical) Turing machines" . 

This idea of classical Turing machines can be enlarged 
to cover Quantum Turing machines, in which the notion 
of a state is generalised to include the usual quantum 
superpositions of states, in the head, tape, register and 
program[4, 5, 6J|. This generalisation seems straightfor- 
ward - however the way one describes such a system has 
led to much misunderstanding and a severe divergence in 
the literature. Computer scientists like to talk in terms of 
states, and logical operations on them - in quantum com- 
puting these operations are called " unitaries" . A " quan- 




FIG. 2: In a quantum walk, the system "hops" between nodes, 
each of which represents a state in the Hilbert space of the 
QUIPS. Allowed transitions are shown as links between the 
nodes. 



turn computation^^ is then a sequence of unitary trans- 
formations Un — Yik ^k performed on some initial state 
Y^in >, to produce the final state |^/ >. The initial 
state is then the input to the computer, and the final 
state the output. 

These states are a quantum mechanical generalisation 
of the binary strings mentioned above - however they are 
much more complicated. Consider, e.g., a classical pair of 
binaries, which can be in one of four states (ie, (11), (10), 
(01), or (00)). In quantum mechanics, systems that have 
2 states are called 'two-level systems' (TLS), or 'qubits'. 
Now however a single qubit can exist in an arbitrary su- 
perposition of states l^i) = ooe^^olo) -|- aie^^^ll), where 
the coefficients ao,ai are real, a^ -\- a\ = 1, the states 
|0) and |1) correspond to the original classical states (0) 
and (1), and </)o, 4>i are phases. A pair of qubits has the 
general wave-function 
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and we see that if we have N qubits then we will end up 
with 2^ coefficients and 2^ 'relative phases' (of which 
one is redundant). The great power of quantum compu- 
tation resides in all the information that can be stored in 
these phases, allowing one to 'superpose' many compu- 
tations at onceQ. A key feature of these states is that 
typically the TV qubits are ^entangled'; one cannot write 
them as a product over independent wave- functions for 
each individual qubit. As first noted by Einstein [3], this 
leads to very counter-intuitive features - it is simply not 
correct to assign a quantum state to an individual qubit, 
and only the A^-qubit wave-function, existing with all its 
relative phases in a vast 2^-dimensional Hilbert space, is 
physically meaningful. 

The above description is very abstract. To get a more 
intuitive understanding, it is helpful to map the QIPS 
onto a "quantum walk" system, in which a single particle 
hops around a graph. Each node of the graph represents 
a state in the Hilbert space of the QIPS, and the links 



between the nodes represent the transition matrix ele- 
ments (in general time-dependent) between them. Thus 
the time evolution of a QIPS can be represented as a 
"quantum walk" in information space [8]; and it is always 
possible to make this mapping for a QIPS composed of 
qubits. 

To see how this works one needs to specify a Hamilto- 
nian, to give the system dynamics. The simplest quan- 
tum walk Hamiltonians take the form[9|: 

Hqw - J2 [Tf^'^iH^- + ^-c-] + E ^M(i)c|.c^ (2) 

where the operators cj^, c^, create or destroy a parti- 
cle at the /x-th node, e^t(i) is the /i-th node energy, and 
Tp,i,{t) the non-diagonal hopping matrix element between 
nodes /i and ly. This Hamiltonian must be generalized to 
describe "composite quantum walks" (where discrete in- 
ternal "spin" variables at each node can also couple to 
the walker) if we are to fully describe a QIPSQ- The 
great advantage of this representation, at least for physi- 
cists and chemists, is that we have a good intuitive un- 
derstanding of how particles move around such graphs. 
One also clearly sees where things can get complicated, 
if the e^ and T^^ are strongly disordered or if their time- 
dependence is non-trivial. How this works in practice can 
be found in the literature [a, 01- 




FIG. 3: Simplified ball-and-stick representation of 

[Fe802(OH)i2(tacn)6]*^ single molecule magnet, show- 
ing the octahedral iron ions, the oxo and hydroxo bridges 
and the 1,4,7-triazacyclononane (tacn) ligands completing 
the coordination spheres. 



"Peg" molecule [l5|, shown in FiglU in which a core of 
eight Fe'^+, spin-5/2 ions forms a total spin-10 system. 

The simplest (and very commonly employed) model of 
a system like Peg starts by reducing the spin Hamiltonian 
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B. Physical Description of a QIPS 

We now turn to a description of a QIPS of the kind 
favoured by physicists and chemists, one which will al- 
low them to design and make one, written in terms of 
the actual physical components of the system. So let 
us now start from a much more physical model, using a 
Hamiltonian written in terms of the qubits themselves: 



ijQB-E4w-r,+Ef^;^'w^r^i 



(3) 



Here the Pauli operator fj operates on the j-th qubit, and 
both the local fields Bj{t) and the interaction U" (t) can 
be time-dependent. 

Now qubits are real physical objects, and many dif- 
ferent physical systems have been proposed to make 
them. Solid-state examples which have achieved 
some experimental success include spins in semicon- 
ductors and quantum dots[lO|, various designs based 
on superconductors [ll|, vacancy centers in diamond [ij], 
single molecule magnets (SMMs)[ij|, and rare earth 
spins fliil. To focus the discussion here, let us consider 
an example. 

An SMM Example: A "Single Molecule Magnet" 
(SMM) typically has a magnetic core of transition metal 
(TM) ions, surrounded by a protective cage of organic 
ligands; the inter-ion couplings are usually antiferromag- 
netic. A useful (because well-understood) example is the 



written in terms of the individual TM spins (for sim- 
plicity we ignore here terms like Dzyaloshinski-Moriya 
interactions) to a single "giant spin" model of form 



H, 



GS 



K2 SaSp 



A^ Oa^lS^j^S 



(5) 



where we assume all low-energy states have the same to- 
tal spin S = '^j Sj. It is very common in the literature 
to drop many of the higher terms in this - for example, 
in the case of Peg one usually writes the simple quadratic 
biaxial form 

n{S) = -DSl + ESl + Ki{SX + St) - 7S • Hi, (6) 

which gives a '2-weir crystal field potential acting on 
S. Here D/kji = 0.23 K, S/fce = 0.094 K, Ki/kji = 
—3.28 X 10~^ K, 7 = gefiBlJ-o, and g^ ^ 2 is the isotropic 
g-factor of the spin-10 moment. All the higher terms 
~ 0{S^) in ([5]) are thus dropped. This effective Hamil- 
tonian for Peg is used below roughly 5 — lOK; above this 
one certainly has to worry about other states, contained 
in the more general 10-spin Hamiltonian ([?]). 

At still lower energies one can truncate (O to its lowest 
singlet or doublet - in the latter case we have a qubit 
Hamiltonian of form 



Hqb = B{Ho) 



(7) 



where the effective field B{Ho) depends both in the crys- 
tal field potential and the applied field Hq. A very im- 
portant special case of this arises for the 2-well potential. 
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FIG. 4: Variation of the effective transverse field Ao by ap- 
plication a real transverse field H^ at different angles in the 
Feg system. 
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FIG. 5: The states of a 3-qubit system, with 8 states, shown 
at left in quantum walk representation, and at right in the 
qubit representation. 



connection to the quantum Turing machine language. 

We have seen above how to describe a QIPS in terms of 
interacting qubits, of a quantum walker, and of a quan- 
tum Turing machine (and there are other descriptions - 
for example, we have for reasons of space ignored adia- 
batic quantum computation entirely). But how well do 
such toy models match the real world in which physicists 
and chemists work? We now address this question. 



C. Realistic Models of QIPS 



where there is a 'natural basis' for the states, in which 
I t) represents the state localised in one well and | |) the 
state localised in the other. One then often writes 
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(8) 
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so that Bx{Ho) — Ao{Hq); here iJg , Hq are the appHed 
fields perpendicular and parallel to the Feg easy axis. 
We can think of Ao as an "effective transverse field", 
driving transitions between the | t) and | J,) states of 
the qubit. This effective field can be varied by applying 
the real external transverse field H{^ - how it does this is 
shown for Fes in FiglH The 2-well system possesses a real 
advantage - it is protected against external transverse 
field fluctuations, which have a small effect on Aq. 

Connection between Quantum Walk and Qubit 
formulations: We have seen that a c^uantum computer 
can be thought of as a particle 'walking' in information 
space, described by a Hamiltonian like Q, or as a set 
of interaction qubits with Hamiltonian ([3]). So how do 
we connect these two? The answer is shown in Fig. [5] for 
a set of 3 interacting qubits - in general the 2^ states 
of an A^-qubit system map to the corners sites of an A^- 
dimcnsional hypercube, and the quantum walker hops 
from one site to another on this hypercube. 

How it hops depends on the exact details of the qubit 
Hamiltonian. Although it is too lengthy to give here, 
the explicit mapping between the quantum walk Hamil- 
tonian Hqb in ^ and the qubit Hamiltonian Hqw in 
([2]), for general qubit couplings, is fairly straightforward - 
it was given in ref. [9| . Another commonly used represen- 
tation is in terms of " quantum gates" [6| , which can also 
be mapped to quantum walks j9,], and thus also allows a 



Obviously both the general models ^, ^, as well 
as the specific models Q, ([S]), and ([5]) for magnetic 
molecules, are oversimplified. Many papers have ap- 
peared discussing both the failures of general models like 
, ([3]) , or of the simple qubit forms ([7]) and ([5]) , as well 
as more specific corrections to the giant spin model ([5|). 
Let us begin by enumerating some of the problems of such 
models as applied to SMM's, and then discuss the more 
general case. The principal shortcomings are as follows: 

1. The giant spin model does not include many other 
discrete electronic degrees of freedom (spin and 
charge); if these lie at low enough energies, they 
couple to T (and of course to S). Because they are 
discrete, they behave like "loose spins" coupHng to 
r (see refs.[la, [l3)- This happens, for example, 
if one has structural transformations in the SMM 
(e.g. Jahn- Teller distortions), or if some exchange 
couplings Jij are much weaker than others. Note 
that in a multi-SMM system, there will also be im- 
purities, and "rogue molecules" which are different 
from the usual SMM because of some internal dis- 
torsion, or because of stray local fields created by, 
e.g., dislocations. These will also behave like "loose 
spin" degrees of freedom coupling to f. 

2. In reality t will couple to many other 'environmen- 
tal modes' in the vicinity of the SMM. This includes 
any nuclear spins in the SMM, coupled via hyper- 
fine interactions, nuclear spins outside the SMM, 
coupled via magnetic dipolar or possibly transfer 
hyperfine interactions [1^, [l3|. Also f will cou- 
ple to photons (via magnetic dipole and electric 
dipole couplings) and both extended phonon modes 



(i.e., magnetoacoustic couplings pTl [JSj ) and local 
phonons (e.g., vibrons, or impurity-induced local 
phonons). In a conducting system, r can also cou- 
ple to extended electronic excitations in various 
ways [131 • 

3. Finally, an ensemble of SMMs, as in a QIPS, will 
have various 'stray' couplings between the {tj}, 
mediated by all the degrees of freedom discussed 
above, and not included in the 'bare' interactions 
U^f^{t) which already appear in ^ above. These 
include 'slow' couplings like those mediated by 
phonons, in which retardation effects can be im- 
portant. 

Actually there is nothing specific to SMM's about 
problems like these. Quite generally, in describing a 
QIPS, we need to include (i) coupling of the qubits to 
stray environmental degrees of freedom, and (ii) stray 
couplings between the qubits from these environmental 
degrees of freedom. 

So what kind of model should we really be using for a 
set of coupled SMM qubits, instead of ([3])? The answer 
proceeds in 2 stages. 

(i) Environmental Modes: We can divide all 
the environmental modes into 2 classes [l^- First, for 
all the discrete localized modes, like defects, nuclear 
spins, "loose" spins, localized phonons, etc, we can use 
a model which couples the qubits to a "spin bath" 
environment [19] of localized modes. In its simplest form 
this describes a set {&%} of two-level systems, and we can 
generalize the bare qubit Hamiltonian in ([3|) to 

Hqb = H^QBir) + Hlsiok) + Hi^in, {ak}) (9) 

in which the Pauli matrices {ffk} act on the 2-level sys- 
tems, and where the spin bath has its own Hamiltonian 



HhiWk}) ^Y.^h,- a, + Y, V^k^lai 



(10) 



fcfc' 



in which the local fields {/ifc} typically dominate over the 
rather weak interspin coupling V^l, ; and finally the qubit 
couples to the spin bath via the term 



fffJ(n,R}) = E-fc"'-"- 



(11) 



with couplings tj^ . More generally we have to couple 
to a spin bath in which the individual modes have more 
than 2 levels - for example, to a set of nuclear spins with 
spin J, having 2/ -I- 1 levels each. This description may 
look rather complicated, but note that we can also write 
Hqb in the more transparent form 

iiQB = H^Qsir) + E (^^ + 6) • ^fe (12) 





qubit network 



FIG. 6: A network of qubits interacts via interactions Jy{t). 
However each qubit is also interactiong with a "spin bath" 
and an "oscillator bath" (for details see text). 



where we define the time-dependent vectors 7fc(t) and 
S,k{t), having components 



kit) = ht + Y^Lufr^it) 
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k' n 



(13) 



Equation (|T2|) simply says that the fc-th bath spin moves 
in a dynamic field coming from both the qubit (incorpo- 
rated in 7fc(t)) and the other bath spins (incorporated 

in ^fc(t)). Note that 7fc(t) is the sum of a static field hk 
and a dynamic field coming from the qubit. Of course 
the motion of the fc-th bath spin (and all the other bath 
spins) will also react back on the qubit via the interaction 

Pa 
^k ■ 

As an example consider the Fcg molecule again, in 
which the core of eight Fe'^"'" spin-5/2 ions couples to 
over 200 nuclear spins in each molecule (these being 
a mixture of ^H, ^H, and possibly ^^O and ^^Fe spin 
1/2 nuclear spins, as well as some higher spin Br and 
N nuclear spins). In this system, analysed in detail in 
refs[20, |2l|, the {V^^,} are weak internuclear couplings 
(l^fcfe' I < 10~^iir even for nearby nuclei), the {w^} are 
hyperfine couplings (with \ujk\ —i-^mK for ^^Fe nuclei 
or nearby proton spins) and the {hk} include both the 
weak nuclear Zeeman couplings, and a term arising when 
the Feg giant spin is not flipping between antiparallel ori- 
entations when it tunnels. 

What about delocalized modes, like phonons or pho- 
tons? These can be handled in a standard way, by cou- 
pling the qubit to a bath of oscillators [22| , which now 
adds an oscillator term 



'Jo.(W)=ii:(^ 



+ mgU^gXl (14) 



plus a coupling 



Hgf{T,{x,})^Y.''>" 



(15) 



to the qubit. In a SMM system, Xq would be e.g. the 
coordinate of a phonon with momentum q, Uq its energy. 



and c" is related to the magnetoacoustic coupling [23J be- 



tween the SMM and this phonon. 
the form 



This coupling takes 
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(16) 



where Os{S) is an operator acting on the giant spin 5, 
0^{xq) an operator acting on the phonons, and g^ a cou- 
pling function; typically l^g-] -- 10 - 100 K for a SMM. 
When we then truncate the giant spin Hilbert space down 
to the qubit space, we get a coupling of the form in (fT5|) . 

Stray Interactions: Now consider the interaction 
U" (t) between the qubits. Essentially this incorpo- 
rates all interactions mediated by fields which operate 
so fast that they can be treated as instantaneous on 
the timescale of the QIPS. This includes all photon- 
mediated interactions (including dipolar interactions and 
exchange interactions). However it does not include 
phonon-mediated interactions, which are slow - if the 
time tc ^ Cs/L taken for phonons to travel across a QIP 
of size L at velocity Cs is not much less than the timescale 
of switching or operation of qubits, then we must treat 
these interactions as retarded, and this can actually lead 
to decoherence. We ignore this problem in what follows 
but it is important. This spin- phonon interaction typi- 
cally has a complicated form - however the main part will 
be of dipolar form, of strength ^^ g'^/pc^rfi, where go is 
the g = component of g^ above. 

It is essential to distinguish between the short-range 
and long-range parts of U^{t); we write 



Utf{t) = K\ 



Q/3 



d: 



Q/3 



(17) 



where Kij is the short-range interaction (typically ex- 
change or superexchange) , with strength ^ 50 — 1000 K 
for transition metal spins and ~ 0.1 — 5 K for rare earths; 
and Dij is the magnetic dipolar interaction, of strength 
^ Vo/rfj, where Vd ~ (S'oo)^ K, where Gq ~ lA is a 
microscopic lattice length. The key point is that Kij is 
only appreciable for nearest-neighbour systems - it falls 
off exponentially and very fast with distance. 




FIG. 7: The k-th bath spin coupled to a qubit has a field 
7fc(t) acting on it, a sum of a static field hk, and a part which 
varies rapidly from — a4 to uJk when 7fe(t) moves between 7^ 
and 7^; this happens when the qubit flips from | i> to | t> 
(a typical path for 7fe(t) during this flip is shown). 



oscillator bath models, but at low temperatures most de- 
coherence comes from spin bath environments of localised 
modes. 

Our purpose in this section is (i) to give readers an 
intuitive feeling for how decoherence works in practise in 
a system of qubits, and (i) give some results for deco- 
herence rates which we hope will help to better design a 
QIPS. 

Quite generally, decoherence in the dynamics of 
some system A is caused by its entanglement with its 
environment |24{ . A key result of quantum mechanics says 
that if we do not know what the environment is doing, 
then we must average over all of its possible states when 
calculating the dynamics of A. But if the environment is 
entangled with A, then this will "smear out" the phase 
dynamics of A. Another way of putting this is that if the 
the environment is reacting to the dynamics of A, then 
it is actually in effect measuring it - and this actually 
destroys phase coherence in this dynamics. 

To see how decoherence works, it is best to divide the 
discussion between the 3 main mechanisms. 



III. DECOHERENCE 

Decoherence is the fundamental problem blocking the 
manufacture of any QIPS. It is considered by many to be 
an insuperable problem (an opinion we do not share!), as 
well as being at the very heart of our understanding of 
quantum mechanics [24]. 

For anyone who wants to make a QIPS, decoherence 
creates a very practical problem - since a QIPS will not 
work if there is too much decoherence, then in design- 
ing one, we need to know how big the decoherence will 
be. This question has a chequered history. Early calcu- 
lations of decoherence were far too optimistic - they used 



A. Spin bath decoherence 

To understand how spin bath decoherence works we 
need first to see how the bath spins move in the presence 
of the qubit. Recall that the k-th bath spin moves in a 
dynamic field 7fe(i) +^k{t) (see equations (IT2l) and (fT3)) '). 
Since ^k is very small compared to 7^, let us begin with 
7fc, the net field on the k-th bath spin (7^; recall it is 
the sum of the static field hk and the time- varying field 
from the qubit. Now 7^ spends almost all of its time in 
the quiescent states 7^. and 7^:, associated with the qubit 
states I t> and | i>, and given by 



7fe = hk + (TLOk 



(18) 




FIG. 8: A spin qubit flipping between its up and down 
quantum states also flips the field acting on nearby nuclear 
spins between two orientations (shown as arrows). The nu- 
clear spins try to precess in this qubit field, but each time 
it suddenly changes they must begin anew. Thus, the path 
a nuclear spin vector follows (an example is shown in red) is 
conditional on the specific trajectory of the qubit, ie., the two 
are quantum-mechanically entangled. 



where a = ± corresponds to | t>, | |> and the vector 



uJk is defined by its components w^ 



(see Fig. [T]). 



However every time the qubit flips (taking a very short 
"bounce time" tb = l/^o to do so), the vector 7fc moves 
rapidly between these two quiescent states. 

The vector ^^ defining the held from the other bath 
spins also varies in time, but quite differently. Typically 
^k will look much like a small and slowly- varying random 
noise source, adding a small "jittering" perturbation on 
7fe, to give the total time-dependent field acting on ak- 

How does CTfc itself move in response to this, as jk 
"jerks" back and forth between 7^ and 7^? The answer 
is shown in Fig [51 At t = 0, ak will begin to precess 
about the initial field ^kit = 0) (which is, say, 7]^). How- 
ever as soon as 7^ jumps to 7^, then ak{t) must begin to 
precess around the new field. The resulting motion of ak 
depends on the exact time sequence of flips of 7^. 

Actually there are two small corrections to this picture. 
First, there is a small extra "wobble" caused by £,k{t)- 
Second, we have described the motion of 7^ as though the 
jumps between 7^ and % were instantaneous, so that ak 
cannot follow them. However this is not true - although 
they are sudden, there is still a small perturbation on the 
dynamics of ak- 

Now notice a crucial consequence of this picture - it is 
that the path followed by ak is conditional on that fol- 
lowed by 7fc . This means that the 2 systems are entangled 
- indeed rather strongly so, since the path followed by ak 



depends more or less entirely on the time evolution of 7^- 
(apart from the small perturbation coming from ^k)- 

Now we can see how a spin bath will cause decoherence 
in the dynamics of a qubit. If we cannot follow the cou- 
pled dynamics of the qubit and each individual bath spin, 
we must average over the different possible time evolu- 
tions of them - a deflnite time dynamics then gets con- 
verted to a time-evolving probability distribution. Thus 
this entanglement between qubit and spin bath dynamics 
causes decoherence in the qubit dynamics. In fact we can 
classify three different decoherence mechanisms here, as 
follows: 

(i) Precessional Decoherence : The precessional motion 
of ak around the time- varying jk , shown in Fig [51 gives 
the strongest decoherence from a spin bath. Essentially 
the spin bath is absorbing phase from the qubit, and it 
turns out that the rate at which the qubit phase dynam- 
ics is 'smeared' is roughly proportional to the average 
solid angle swept out by all of the ak together in their 
precessional motion. A detailed calculation [l^ [20| gives 
simple results for the decoherence rate F^ in two limiting 
cases: 

F^ -- l/2^(wfe//ife)2; (iUk^hk) (19) 

fc 

Fj ^ l/2Y,ihk/^kf; (hk-^tok) (20) 

k 

where cuk = \uJk\ in (|18p . and hk = \hk\- Note that this 
"precessional decoherence" (originally called "orthogo- 
nality blocking" [25| ) does not involve dissipation -no en- 
ergy is exchanged between qubit and bath, only preces- 
sional phase. Thus this decoherence mechanism, by far 
the most important one coming from the spin bath, is 
not even dissipative! 

The simplest and best-studied example of precessional 
decoherence is the case of a spin qubit coupled to nuclear 
spins. Consider the example of Fcg, already discussed 
above, where the main contribution to hk is just the Zee- 
man coupling to the external field, i.e., hk ~ gNfJ^NlkHo, 
where gpf and /zat are the nuclear g-factor and magnetic 
moment. Then from (|20)) we see that when Hq = 0, the 



decoherence rate F^ = (this is because 7^ flips then 
through 180°; such a flip has no effect on the dynamics 
of (T/c). Decoherence then rises rapidly until reaching a 
maximum when hk — uJk, where cuk is the hyperfine cou- 
pling as before - at this point the bath spins are precessing 
wildly as the field jk flips between orientations separated 
by angles ~ 90°. Further increase in Hq then decreases 
F^, and in the high-field limit F^ again becomes very 
small (in this case ^l and 7^ are almost parallel, so their 
flipping of Tz again hardly affects the dynamics). The 
quantitative details for Fes were worked out in ref . |20j . 

(ii) Topological Decoherence: The above discussion 
concerned the bath spin dynamics when the qubit was 
sitting quietly in either | t) or | J,). But what about 

when the qubit is actually flipping, and 7fc(i) is moving 

I —ft 

rapidly between 7^ and 7^? If this flip were instanta- 



neous, there would be no effect on the bath spins - they 
simply couldn't follow the sudden change. However the 
time 1/rio taken for jk to flip from 7^ to 7^ is actually 
finite (though still very short). Then elementary time- 
dependent perturbation theory tells us that the ampli- 
tude ak for the spin ak to make an inelastic transition 
under this sudden change, ie., during the qubit flip it- 
self, is ak ~ 7rcjfc/2r2o- By the same arguments as before 
this bath transition must cause decoherence in the qubit 
dynamics - the bath is actually registering the change in 
the qubit state (it is in effect performing a measurement), 
and this transition actually adds a random topological 
phase to the qubit dynamics [lg|. The resulting "topo- 
logical decoherence" contribution FT to the decoherence 
rate is found to be[la.[l9|: 



1/2 E 



ak\ 



(21) 



M=0 contribution 



0.2 - 



0.1 



For most systems this will be very small, since ujk ^ ^0 
as a rule. Thus, for Feg, fio ~ 5K in low applied fields, 
whereas most of the lo^. < ImK; even after summing over 
c. 200 nuclear spins, one still finds a very small contri- 
bution to decoherence, changing slowly with field as flo 
decreases f2(|. 

(iii) Noise Decoherence : The best-known kind of deco- 
herence is that coming from some unknown noisy vector 
field x{t) coupling to r. Chemists are very familiar with 
this in NMR and EPR - it leads to the usual "T1/T2 phe- 
nomenology" , in which the noise causes incoherent relax- 
ation of the qubit energy over a time Ti and the qubit 
phase over a timescale T2. A spin bath also gives rise to 
this 'noise decoherence'; the weak interactions between 
the bath spins cause a slowly fluctuating field to act on 
the qubit, causing Ti and T2 relaxation in the standard 
way. 

Notice an interesting consequence of this, which is seen 
very nicely in experiments on SMMs. Since most of the 
spin bath dynamics is being driven by the qubit itself 
(see again Fig |5]) , then when the qubit is frozen by an 
external longitudinal bias field e ::^ A, the bath spins 
must also largely freeze - without the time- varying field 
7(i) to drive them, only the very weak Vkk' are left to 
give the spin bath any independent dynamics (via the 
small slowly- varying (.kit))- Thus we expect Ti and T2 
to drastically increase when the qubit is frozen. This 
was seen in very pretty experiments by the Leiden group 
on the SMM Mni2; they observed roughly a hundredfold 
increase in Ti and T2 when they pushed the tunneling 
Mni2 molecules off resonance f20|. 

Further remarks : An obvious questions which arises 
here is the following: to what extent are all these deco- 
herence results already incorporated into the T1/T2 phe- 
nomenology used in experimental chemistry and physics 
in the last few decades? 

We note first that the common practise in the literature 
(particularly the quantum information literature) of asso- 
ciating decoherence with some "experimental noise" fluc- 
tuation source, and then applying standard ideas such as 
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FIG. 9: Plot of Im xtt('^)i the absorption probability de- 
scribed in the text, for a qubit coupled to a spin bath. With- 
out decoherence this would be a sharp line at lj/Ao — 1. 
The shaded part shows contributions from processes where 
the spin bath polarisation is unchanged. 



the fluctuation-dissipation theorem[27| to associate this 
noise with a dissipation, is quite wrong for a spin bath. 
As noted above, the most important part of spin bath de- 
coherence, i.e., precessional decoherence, simply involves 
precessional motion - there is no transfer of energy be- 
tween qubit and bath spins, and therefore no dissipation. 
Thus there is no classical analogue to decoherence from 
a spin bath, and no connection between decoherence and 
dissipationlll [H, [H [H . 

Now the "T1/T2" phenomenology in NMR is based on 
the classical Bloch equations. It throws away all the in- 
ternal phase information in the nuclear spin bath dy- 
namics (including relative phases between spins, entan- 
glement, etc), leaving only the self-consistent description 
of a single nucleus in the classical fluctuating mean fleld 
from all the other spins (plus any fluctuating fields com- 
ing from electron spins, fluctuating external flelds, mo- 
tional narrowing, etc). 

However this T1/T2 phenomenology, although it is 
commonly used in the qubit literature (in discussions of, 
e.g., Rabi oscillation, or spin echo for a qubit), it is not 
appropriate to the qubit dynamics [l9|, |28|. To give an 
idea of what can really happen, in Fig. [9] we show the 
imaginary part of the Fourier transform X]]{^) of the 
spin correlation function P]-\{t) for a qubit coupled to a 
spin bath. This function is the probability that a qubit 
in state | t> at the time i = will be found there at time 
t. Thus X||(a;) is the absorption intensity at frequency 
w in a qubit resonance experiment. This looks nothing 
like the Lorentzian expected in the T1/T2 phenomenol- 



ogy, and it shows singular peaks which come from the 
unusual coupled dynamics of qubit and spin bath. These 
features actually result from decoherence and are not as- 
sociated with dissipation. 

This point is important, because it means that eventu- 
ally experimentalists working on qubits are going to have 
to abandon the T1/T2 phenomenology in describing their 
experiments, when spin bath decoherence is important. 
The phenomenology that will be required needs a whole 
paper on its own. 



B. Oscillator Bath Decoherence 

This has been well covered in other articles [2^ [23| and 
books [30]. The theoretical problem is simply formulated 
- one now couples a qubit to a bath of oscillators, and 
tries to calculate th e dy namics. This "spin-boson" model 
is well understood|3(3|. In contrast to the spin bath, 
there is a quantum-classical correspondence for oscilla- 
tor baths [3l|, [34I, and a fluctuation-dissipation theorem 
relating both the dissipative dynamics and the decoher- 
ence of the quantum system to the fluctuational " quan- 
tum noise" acting on it, coming from the oscillator bath. 

Thus we can understand decoherence in the spin-boson 
model in a fairly simple way. The oscillator bath has two 
main effects. First, it slows down the qubit, increasing its 
"inertia", so that A, the qubit flip rate, is " renormalizcd" 
(ie., A ^ A with A < A). Second, it causes dissipation 
in the "quantum relaxation" of the qubit dynamics, ei- 
ther because the qubit spontaneously emits an excitation 
(i.e., excites an oscillator from its ground state) or scat- 
ters an existing excitation (i.e., stimulated emission - an 
excited oscillator changes its state, and possibly another 
one is excited). All the complexity in the qubit dynamics 
comes from repeated absorption and emission of oscilla- 
tor excitations |23. ISSj . Nevertheless the phenomenology 
is very familiar - one defines a "spectral function" '31] of 
form 



Ja{Lu)^n/2j2- 



r.a|2 



-5{UJ - LUq) 



(22) 



i.qujq 



(where a =||,_L), and the entire dynamics can be cal- 
culated in terms of this function, which is familiar from 
second order perturbation theory or from Fermi's golden 
rule. It has the form (I/^/w)x (density of states), where 
V^ = c" is a coupling constant and m the energy denomi- 
nator, and the sum over 5{lo — coq) measures the density 
of oscillator states at uj. 

Most of the theoretical work on the spin-boson system 
has been done for the simple T- independent "Ohmic" 
form J{lu,T) — j rjcu. However most experiments, cer- 
tainly on qubit systems, are not described by an Ohmic 
model, for two reasons. First, as already discussed, spin 
bath decoherence is not at all described by such a model. 
Second, even without the spin bath, the oscillator bath 
is often not Ohmic. It certainly will not be if the system 



is insulating, since for phonons J(w, T) depends strongly 
on T and is a higher than linear power of uj. Only for 
simple conductors is the Ohmic form always valid. For 
a proper discussion of the different forms of J{u!), see 

As an example let us again take the SMM. The prin- 
cipal source of oscillator bath decoherence is then the 
magnetoacoustic coupling to acoustic phonons (for which 
J{uj) ~ cj'^. Both quantum relaxation rates [13, [3 and 
decoherence rates [20, [SJ] have been calculated for this 
interaction. The usual form for the decoherence rate in 
a transverse field is 



\MffA^^ . .A 



nn^pcf''''^^kf^ 



(23) 



where M/. is the spin-phonon matrix element discussed 
in section 2, p the density, and Cs the sound velocity. 
This expression illustrates nicely the points made above; 
we have a (matrix element)^, multiplied by a density of 
states ~ Aq, divided by an energy denominator Ao, mul- 
tiplied by a thermal factor ~ A"-'^ coth(Ao/fcr). In other 
words, just the golden rule. A more sophisticated calcu- 
lation hardly changes this result, simply renormalizing 
Aq. An Ohmic bath would give a completely different 
result [13, [221 • Actually, decoherence from an Ohmic bath 
is typically very large, which is why in making a QIPS we 
want to avoid having moving electrons (ie., one should 
use insulating systems). Of course, this still leaves de- 
cohrence from the spin bath. 



C. Pairwise Dipolar Decoherence 

We now turn to what may in the long run be the most 
insidious kind of decoherence - that caused by unwanted 
long-range interactions between the qubits. There are 
two kinds of these, viz. (i) interactions like dipolar or (in 
the case of current-carrying devices) unscreened induc- 
tive interactions, which can be treated as instantaneous; 
and (ii) long-range interactions mediated by slow parti- 
cles like phonons, which may have a retarded character. 

The reason such interactions are so dangerous is that 
they affect many qubits at a time. Only a few theoretical 
studies have so far appeared. Aharonov et al.[33 argue 
that concatenated error correction codes can handle long- 
range interactions, but only just. It is however likely that 
the model in this paper is over-simplified, since detailed 
calculations for both insulating systems [34l| (where dipo- 
lar inter-qubit interactions are important) and conduct- 
ing systems [36| (where inter-qubit interactions mediated 
by the electron bath are important) give rather different 
results. Thus the problem is still open, and we do not 
give any detailed results for decoherence rates here. 

We can nevertheless give a simplified description of 
what is going on (see Fig. [TUl). With long-range (~ 1/?'"^) 
interactions, it is always possible for a given qubit to find 
resonance with many others in 3 dimensions, whether 
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FIG. 10: 2-dimensional spin network, showing resonant pairs 
of qubits interacting via dipolar couplings. 



we want it or not, provided the strength of these inter- 
actions exceeds a critical threshold (this result basically 
goes back to Anderson|37t). Typically for dipolar inter- 
actions the resonances will occur over rather long dis- 
tances. There are really only two ways to get round this 
problem - first, to use error correction (a topic whose 
details [6| we have eschewed in this article); and (ii) use 
a lower dimensional geometry, which drastically reduces 
the number of qubits that are far apart (compared to the 
3d case). Ideally one uses both - in section 4 we dis- 
cuss certain architectures which can reduce the effective 
dimensionality of the QIPS. 



IV. DESIGN CRITERIA: POSSIBLE 
ARCHITECTURES 

Let us now use the results given in the previous section 
to discuss the best way to make a QIPS. To recapitulate 
- we want a system of qubits whose mutual interactions, 
and the local fields acting on them, can be accurately con- 
trolled (or, if they are static in time, are very accurately 
fixed and can be tuned). We want to be able to write 
quantum information onto the states of these qubits, do 
a computation, and then read out the answer. We do 
not want stray interactions or fields interfering with this, 
and we certainly want to minimize decoherence caused by 
interactions with other uncontrolled degrees of freedom. 
So how can we best achieve all these desiderata? 

As advertised in the introduction, we concentrate here 
on designs for spin QIPS using materials chemistry, ie., 
'bottom-up' approaches. The main advantages of such 
an approach (as opposed to top-down 'nanofabrication' 
approaches), are: 

Al) reproducibility - chemistry uses quantum mechan- 



ics to produce absolutely identical molecules or other 
kinds of spin complex; 

A2) Size - the qubits and other parts of the architec- 
ture are really small, so that few degrees of freedom are 
involved, and few defects; hence decoherence is reduced 
and timescales become short; 

A3) No moving electrons - moving electrons always 
cause decoherence. With insulating molecules or other 
spin complexes one can make a true spintronics system, 
which only involves spin degrees of freedom. 

There are also disadvantages: the main ones are 

Dl) Size - molecular and ionic spin moments are cur- 
rently too small to be either controlled or examined in- 
dividually, except in very small numbers. In most QIP 
algorithms, we need to be able to control, address, and 
read out the states of individual qubits. 

D2) Dipolar interactions - as discussed above these are 
potentially fatal to QIPS - one needs architectures that 
suppress them. 

Our task is thus to find designs that use the advan- 
tages (A1)-(A3) to their fullest extent, and overcome the 
disadvantages (Dl) and (D2) as much as possible. In 
bottom-up designs the main weapons we have are (i) the 
methods of chemistry, and (ii) a huge variety of potential 
geometries, ie., possible QIPS 'architectures'. We deal 
with these in turn. 



A. Chemical approaches 

At the level of individual ions or molecules, materials 
chemists are well-placed to deal with the problems of 
reproducibility and decoherence, both at the nuclear and 
electronic scale. 

(i) Nuclear Chemistry: Unless one wishes to use 
nuclear spins as a QIPS resource (for which see below), 
then their enormous contribution to decoherence means 
we simply want to get rid of them. Now while all "odd" 
nuclear isotopes (ie., with an odd number of protons) 
have / 7^ 0, at least one isotope of every even element 
has / = 0. Indeed, in natural mixtures of even-numbered 
elements, the / = isotopes are the most common. Thus 
we should aim for systems with even-numbered elements 
only. 

Of course, if one aims strictly for even nuclei, the chem- 
ical playground is dramatically reduced, and isotopic pu- 
rification is still typically required. Obvious transition 
metal choices are Fe (with 2% ^'^Fe) or Ni (with 1% ^^Ni). 
There are also a number of rare earth candidates. These 
are either even-numbered (like Nd, Dy, Er, Yb) with a 
complex mixture of isotopes and the possibility of having 
only / = through costly isotopic purification, or odd- 
numbered (like Tb, Ho), where / y^ but only a single 
isotope is naturally prevalent. Note however that in this 
case we deal not with a single rare earth ionic doublet at 
low energies, but an electronuclear pair of spins - the nu- 
clear spin hyperfine coupling is so strong that one cannot 
rotate without the other I38l. 
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(a) 



FIG. 11: Some systems that can be synthetized with a smaU 
amount of nuclear spins and that could serve as hardware for 
qubits (see text for details): (a) Single-molecule magnet based 
on a single lanthanide ion. (b) Bidimensional "honeycomb" 
oxalate network. 



An entirely new range of possibilities opens up if we 
use nuclear spins as part of the QIPS (e.g., in a quan- 
tum register). We then must be able to switch on and 
off the interaction (and hence the information flow) be- 
tween the nuclear and electronic spins. One way to do 
this is dynamical - one tunes the electronic and nuclear 
spin splittings (e.g., with external fields) so that resonant 
hyperfine-mediated interactions can be switched on and 
off. Concrete practical designs are needed here - the big 
stumbling block will be decoherence from dipolar interac- 
tions between nuclear spins and remote electronic spins 
(even if these latter do not partake in the resonant inter- 
actions, their dipolar field adds extra dynamic phase to 
the nuclear dynamics, and can throw them off resonance). 

Another possibility is to switch the hyperfine interac- 
tions on and off. In spin transition systems the electron 
spins can be switched on and off without a change in 
the redox state, using temperature, pressure or light (the 
LIESST effect, m or HAXIESSTf4li / SOXIESST^ 
for hard / soft X-rays). If a design based on this effect 
could be made to work it would be very powerful - one 
would simply switch off the QIPS when it was not needed, 
storing information in the nuclear register. 

(ii) Electronic Chemistry: Undoubtedly one area 
where chemists can play a big role is in the design of 
interactions (ie., engineering of effective Hamiltonians) , 
and the preparation of uniform & reproducible qubits. 

The first problem then is is the design of individual spin 
qubits. Let us first note that seemingly obvious candi- 
dates such as molecular Kondo s yste ms or electrically- 
controlled quantum dots[i3, IS Ej], while promising, 
have one very serious drawback - they have an intrin- 
sically high decoherence, as their operation is based on 
moving electrons. 

So let us now consider insulating SMMs and rare earth 
ions as qubit candidates. A question often asked by 
molecular chemists is - what kind of crystal field is needed 
for the spin qubit? As noted earlier, designs in which 
there is a 'natural basis' for the qubit, created by a 2- 
well potential, are usually good ones. Many tunneling 



SMMs and rare-earth ion molecular systems now use this 
design - the question at issue is how big to make the spin 
anisotropy energy barrier Eb- One wishes to operate 
the QIPS at energy scales ^ kTc, where kTc ~ Eb/^ttS 
is the crossover energy to the qubit regime (above Tc, 
the spin dynamics proceeds by incoherent thermal acti- 
vation). However at the same time we want to suppress 
decoherence from nuclear spins and phonons. This im- 
poses competing restrictions: to reduce nuclear spin de- 
coherence one wants Aq » {^k}, where the {tUk} are 
the couplings to the spin bath (nuclear spins in this 
case), but to reduce phonon decoherence one wants Aq 
as small as possible compared to typical phonon ener- 
gies (whose scale is determined by the Debye energy 
9o)- Thus there is a 'coherence window' zO], defined 
by 9d ^ ^o ^ {^k}, defining a desirable size for Aq. 
Noting that Aq is certainly less than kTc, we can put 
these various restrictions together as a requirement that 
Od > Eb 127:3 > kT, A^ > {ujk}- Thus we would like 
to have a large magnetic anisotropy and weak hyperfine 
couplings - and a system that can go to very low temper- 
atures. Certainly the most important requirement here 
is weak hyperfine couplings. It also helps if Od is big and 
if the spin-phonon couplings are weak (which will happen 
if the molecule is, eg., nearly spherically symmetric). 

If one uses transition ion SMMs based on even nuclei 
(so as to eliminate nuclear spin decoherence), severe re- 
strictions are placed on the chemical bridging between 
the spins - one has to discard all the halogens and almost 
every organic molecule because of their hydrogen and ni- 
trogen nuclear spins. For coordination and short-range 
bridging, obvious choices are then oxo, tio or oxalate an- 
ions, but there are others, like carbonyl, thiocarbonyl or 
carbide. 

Alternatively one can go for rare earth-based 
molecules. Single lanthanide cations in adequate coor- 
dination environments have shown SMM behaviour |44|: 
and recently, polyoxometalates (POMs) composed com- 
pletely of W, O and Si have shown crystal fields capa- 
ble of giving qubit dynamics [4^. This offers exciting 
new possibilities, as POMs are chemically quite different 
from the usual SMMs (they have high negative charge, 
and are stable against oxidation). Two advantages of 
any rare earth design are (i) the relevant magnetic pa- 
rameters can be calculated theoretically from their struc- 
tures and/or extracted from EPR and optical data; and 
(ii) their strong magnetic signals and hyperfine coupling 
(when using nuclei with / > 0) facilitate their detection 
and characterization. 

Another key requirement is the reduction of disorder to 
very low levels - this means eliminating defects, disloca- 
tions, etc., from the QIPS, and using molecules or other 
spin complexes that only come in one species (ie., not 
using systems that have active structural degrees of free- 
dom, various kinds of isomerism, or free rotation about 
single bonds, or solvents/counterions that can take dif- 
ferent positions which are close in energy to each other). 

Let us turn now to the design of inter-qubit interac- 
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tions. These need to be precisely controlled (whether 
they are static or dynamic) and inter-qubit dipolar in- 
teractions and spin-phonon interactions need to be sup- 
pressed. Since both are oc the qubit spin moment, one 
possibility is to have antiferromagnetically ordered spin 
qubits. The inter-qubit interactions would then vanish 
much more quickly at long range (for an antiferromag- 
netic pair, they are oc a/r'^, where a is the separation 
between the spins). If the spins in each pair are not 
equivalent (eg., in an antiferromagnetic transition metal 
dimer where the two metals are either different or have 
different chemical environments), an external stimulus 
(e.g. an EPR pulse) will differentiate the | }[> and 
I J.t> configurations. This way of suppressing dipolar 
interactions also suppresses the spin-phonon interaction. 
Important work that has already been done in this di- 
rection with dimeric molecular nanomagnets|46l WA. 4ql. 
including molecules where there is a built-in switch 47| . 
We would suggest making analogues to these systems in 
which the dipolar magnetic moments cancel completely 
while a strong quadrupolar moment remains. Alterna- 
tively, heterometallic antiferromagnetic wheels [4^, with 
S = Q, are clearly promising building blocks, and qubit 
designs for these systems have already been explored [50|. 

The most obvious way to implement inter-qubit in- 
teractions is to use exchange or superexchange inter- 
actions between neighbouring spins. Transition metal- 
based systems have interactions ^ 0{1QQK), making 
them better than rare earths, where these interactions 
are typically ~ 0{1K). If controllable long-range inter- 
qubit interactions are desired, one can imagine propa- 
gating these through reduced carbon nanotubesfSl] or 
POMs^]. Both can in principle mediate relatively long- 
range indirect exchange through the delocalized elec- 
trons, compared to superexchange, and both can be pre- 
pared using only / = isotopes. However, they involve 
moving electrons, a serious serious source of decoherence; 
a better way would be to propagate spin signals down a 
chain of strongly exchange-coupled insulating spins. 



B. Geometry, Architecture, and Fabrication 

The actual spatial arrangement of qubits, and the 'read 
in/read out' probes which couple to them, is quite cru- 
cial to the design of a QIPS. There are 3 main issues, 
viz. (i) the small size of spin qubits based on molecules 
or other spin complexes means that with current detec- 
tion/control systems, it is hard to manipulate, read in, or 
read out the state of a single qubit; (ii) the much larger 
size of these detection/control systems means that they 
cannot be 'crowded in' to interact with more than a few 
qubits at a time; and (iii) the geometrical design of a 
QIPS will greatly influence decoherence, both from stray 
fields (e.g., from substrates) and from long-range dipolar 
interactions. 



(i) Geometry and Architecture: Current probes 
(microSQUIDs, STM and MFM systems, HaU probes, 
optical detection and control systems, etc.) are partly 
dogged by the uncertainty principle - if one is not careful, 
attempts to probe or control at a length scale of order 
a qubit size (ie., a few nm) will have far too destructive 
an effect on the qubits and their quantum states (thus 
if photons were used, we would be dealing with soft X- 
rays!). STM and MFM probes can have a very weak 
effect, even though their tips are very small - but the 
rest of them is very big, and there is no obvious way to 
address many qubits at a time with this technology. 

The most obvious way to get round this problem is 
to use fixed probes of nm size, connected remotely to the 
outside world along one-dimensional connecters. Conven- 
tionally one thinks here of wires and moving electrons 
- conceivably some future nanoSQUID array could do 
the job. However it seems better to us to use a gen- 
uinely spintronic arrangement of interacting spins (e.g., 
using strongly-coupled chains of spins, or possibly a set 
of nanotubes) to couple into the qubit array. These can 
themselves be connected to much larger control and read 
in/read out systems outside the QIPS. As yet there are 
no concrete designs in this area - it would be useful to 
have some. 

Another possibility, at least in the near future, is to 
make larger spin-based qubits - a large (e.g., 50 x 50) spin 
array (square, hexagonal, linear) of microscopic spins or 
SMMs, coupled antiferromagnetically. These would in- 
teract via a few nearest neighbour spins. A further ad- 
vantage of this design would be obtained if each qubit had 
no net spin, thereby suppressing dipolar errors. Promis- 
ing advances here are being made with the s ynt hesis and 
study of mesoscopic antiferromagnetic grids [5a|. 

Another question which we find interesting is the way 
in which both the read in/read out and decoherence prob- 
lems might be alleviated by using a 'sparse architecture' 
for the qubits [56|. By this we mean an arrangement of 
qubits on, e.g., a planar substrate, whose total number 
does not increase linearly with the area of the system, 
but more slowly. This allows easier access for probes, 
and also strongly reduces errors caused by dipolar in- 
teractions. Many geometries are possible here, ranging 
from coupled lines or nanorods of various shapes to 'frac- 
tal' patterns on surfaces. The chemistry and physical 
properties of the substrates will be important in limiting 
these designs. For example, while dendrimers[57] could 
be a promising scalable support for sparse architectures, 
steric hindrance usually induces disorder in their struc- 
tures, which would cause decoherence. 

Finally, let us mention two other interesting alterna- 
tive designs which could alleviate the architectural prob- 
lem. The first involves using specific alg orithms that 
avoid single-qubit addressing altogether [53| ; only the en- 
tire QIPS is addressed, but on multiple occasions after 
a series of pulses. The other uses molecular cellular au- 
tomata [5J| which only need addressability at edges, thus 
reducing the addressability problem by one dimension. 
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Both of theses designs suffer because they require a long 
time to carry out a computation - they can only work if 
decoherence times are really long. 

(ii) Fabrication: This brings us back to the problem 
of fabrication of these QIPS arrays. The range of possibil- 
ities is enormous; we mention just a few. As far as SMMs 
are concerned, we gave the main desiderata above - the 
interesting question is what new possibilities would be 
useful to explore. One is POMs, which have a rich chem- 
istry offering different topologies, sizes, and the ability to 
encapsulate magnetic metal ions in arbitrary ligand fields 
and/or host delocahzed electrons if desired [58]. As noted 
above, they can be built entirely using even elements 
(mainly W, Mo and O). Thus, one could explore "giant" 
POM wheels |59l|. currently about 4nm in radius. It may 
be possible to build a wheel which (i) is even bigger, (ii) 
is free from nuclear spins (current designs contain many 
water molecules), and (iii) has the magnetic network we 
want inside. Alternatively, with an array of STMs, one 
could chemically prepare a self-assembled monolayer of 
POMs or fullerenes, and inject a single electron on each 
molecule directly under each tip. The chemistry would 
be quite easy, and avoids nuclear spins, since fullerenes 
and POMs are based on even-numbered elements. The 
main problems would be (i) electron delocalisation in the 
molecules introduces new paths to decoherence, and (ii) 
this scheme would shift part of the difficulty from the 
chemistry to the nanoengineering. 

We also note that a way to prepare antiferromagnetic 
2-D lattices which are poor in nuclear spins and have 
two distinct sites is to use bimetallic honeycomb oxalate 
layers j6(]f: see Fig. [TT] One of the metals will always 
carry a nuclear spin, but the other can be clean if one so 
chooses. They have no water, but do require a nearby 
cationic layer, to compensate their own anionic charge. 
This cationic layer could be Ca^^-based. 

Finally, we emphasize the care required in designing 
substrates. They will cause strong decoherence if there 
are any 'spin bath' defects (dangling bonds/free radi- 
cals, or charge defects, or paramagnetic spins, or nu- 



clear spins) which can couple to the qubits. Solving 
this will not be easy - getting rid of nuclear spins in 
the substrate may be the hardest task of all. POMs, 
among a vast variety of molecules, have been organized 
bidimensionally in different ways, including Langmuir- 
Blodgett films and through covalent modification of the 
surface. Of course, usually these ways make extensive 
use of elements with an odd number of protons, so po- 
tentially the cleanest way would be to use self-assembled 
monolayers (SAMs) on an adequate substrate. Highly 
ordered pyrolytic graphite 61 1 and silicon|62l| have been 
demonstrated to support SAMs of - among other chemi- 
cal species - different POMs, and are all free from nuclear 
spins. Quartz would be a good insulating substrate, but 
the preparation of POM SAMs on quartz usually involves 
organic molecules to provide for a positive charge ^63j, 
which defeats the purpose of choosing a nuclear-spin-free 
substrate. 



V. CONCLUDING REMARKS 

The international effort to make a QIPS has assumed 
very large proportions, but building a QIPS will obvi- 
ously not be easy - the main problem, as we have seen, is 
to suppress decoherence, particularly from nuclear spins 
and dipolar interactions. The use of insulating atomic or 
molecular scale spin qubits offers important advantages 
- reduced decoherence, and with appropriate care, high 
reproducibility. We have explored herein a large number 
of design strategies for building a QIPS. It seems likely 
to us that experimental efforts to gain control of decoher- 
ence and make systems of multiply-entangled spin qubits 
will succeed in the next few years. At this point it should 
be possible to devise realistic architectures for large scale 
quantum information processing, and we anticipate that 
chemical considerations will play a large role in these ef- 
forts. 

This work was supported by NSERC, CIFAR, and 
PITP in Canada, and by MEG in Spain. 



[1] A review of ion trap work is in R Blatt, D Wineland, [10] 

Nature, 2008, 453, 1008. 
[2] A.M. Turing, Proc. London Math. Soc, 1936, 42, 230 
[3] R.P. Feynman, " The Feynman lectures on computatiori' , 

Addison- Wesley (1996) [11] 

[4] R.P. Feynman, Found Phys., 1986, 16, 507 
[5] P. Benioff, Phys. Rev. Lett, 1982, 48, 1581, and J. Stat 

Phys., 1980, 22, 563 
[6] M.A. Nielsen, I. L. Chuang, " Quantum Computation and [12] 

Quantum Information" (CUP, 2000) 
[7] A. Einstein, B Podolsky, A Rosen, Phys. Rev., 1935 47, 

777 [13] 

[8] see J. Kempe, Contemp. Phys., 2003, 44, 307, and refs. 

therein. [14] 

[9] A.P. Hines, P.C.E. Stamp, Phys. Rev A, 2007, 75, 062321 



T. Fujisawa et al.. Nature, 2002 419, 278; T Hayashi 

et al, Phys. Rev. Lett, 2003, 91, 226804; K Ono et al., 

Science, 2002, 297, 1313; J.R. Petta et al., Science, 2005, 

309, 2180; F.H.L. Koppens et al. Nature, 2006, 442, 776 

Y Nakamura et al.. Nature, 1999, 398, 786; C.H. van 

der Wal et al.. Science, 2000, 290, 773; D. Vion et al., 

Science, 2002, 296, 886; M Steffen et al., Science, 2006 

313, 1423 

F Jelzko et al., Phys. Rev. Lett, 2004, 92, 76401; ibid, 

2004, 93, 130501; T Gaebel et al., Nature Phys., 2006, 2, 

408; MVG Dutt et al. Science, 2007, 316, 1312 

S Hill et al, Science, 2003, 302, 1015; S. Bertaina et al., 

Nature, 2008, 253, 203 

J Brooke et al.. Science, 1999, 284, 779; R. Giraud et al., 

Phys. Rev. Lett, 2001, 87, 057203, and ibid, 2003, 91, 



14 



257204; S. Bertaina et al., Nature Nanotech, 2007, 2, 39; 

S Nettutlaet al., Phys. Rev. Lett, 2007, 99, 137601 
[15] D. Gatteschi et al., Science, 1994, 265, 1054 
[16] N.V. Prokof'ev, P.C.E. Stamp, J. Phys CM, 1993, 5, 

L663 
[17] N.V. Prokof'ev, P.C.E. Stamp, J, Low Temp. Phys, 1996, 

104, 143 
[18] F. Hartmann-Boutron et al.. Int. J. Mod. Phys., 1996, 

BIO, 2577 
[19] N.V. Prokof'ev, P.C.E. Stamp, Rep. Prog. Phys., 2000, 

63, 669 
[20] P.C.E. Stamp, I.S. Tupitsyn, Phys. Rev. B , 2004, 69, 

014401 
[21] I.S. Tupitsyn et al., to be published. 
[22] A.J. Leggett et al. Rev. Mod Phys., 1987, 59, 1 
[23] E. du Tremolet de Lacheisserie, " Magnetoostriction: 

Theory and applictions of magnetoelasticity" (CRC press, 

Boca Raton, 1993) 
[24] For introductory discussionos of decoherence, see R. 

Omnes, " The interpretation of Quantum Mechanics" , 

Princeton Univ. Press (1994); and P.C.E. Stamp, Studies 

Hist. Phil. Mod. Phys. 37, 467 (2006) 
[25] N.V. Prokof'ev, P.C.E. Stamp, pp. 347-371 in "Quantum 

Tunneling of Magnetisation: QTM '94 ", ed. L. Gunther, 

B. Barbara (Kluwer, 1995) 
[26] A. Morello et al., Phys. Rev. Lett, 2004, 93, 197202; A. 

Morello, J de Jongh, Phys. Rev., 2007, B76, 184425 
[27] The fluctuation-dissipation theorem is explained in, e.g., 

P.M. Chaikin, T.C. Lubensky, "Principles of Condensed 

Matter Physics", C.U.P. (1995) 
[28] M. Dube, P.C.E. Stamp, Chem. Phys., 2001, 268, 257 
[29] H. Grabert, P. Schramm, G-L Ingold, Phys. Rep., 1988 

168, 165 
[30] U. Weiss, " Quantum Dissipative Systems" , World Scien- 
tific (1999) 
[31] A.O. Caldeira, A.J. Leggett, Ann. Phys., 1983, 149, 374 
[32] A.J. Leggett, Phys. Rev., 1984, B30, 1208 
[33] R.P. Feynman, F.L. Vernon, Ann. Phys., 1963, 24, 118 
[34] A. Morello, P.C.E. Stamp, Phys. Rev. Lett, 2006, 97, 

207206 
[35] D Aharonov, J Preskill, A. Kitaev, Phys. Rev. Lett, 2006, 

96, 050504 
[36] E. Novais, H.U. Baranger, Phys. 

040501; E. Novais, E.R. Muccioh, 

Rev. Lett, 2007, 98, 040501 
[37] P.W. Anderson, Phys. Rev., 1958, 109, 1492 
[38] M. Schechter, P.C.E. Stamp, Phys. Rev. Lett, 2005, 95, 

267208; M. Schechter, P.C.E. Stamp, Phys. Rev B, in 

press (/condmat 0801.2889). 
[39] S. Decurtins, P. Giitlich, C. P. Kohler, H. Spiering, A. 

Hauser, Chem. Phys. Lett 1984, 105, 1-4 
[40] D. CoUison, C. D. Garner, C. M. McGrath, J. F. W. 

Mosselmans, M. D. Roper, J. M. W. Seddon, E. Sinn, N. 

A. Young, J. Chem. Soc. Dalton Trans. 1997, 4371-4376 
[41] G. Vanko, F. Renz, G. Molnar, T. Neisius and S. Karpati, 

Angew. Chem. Int Ed. Eng., 2007, 46, 5306-5309 
[42] L. Gao, W. Ji, Y. B. Hu, Z. H. Cheng, Z. T. Deng, Q. 



Rev. Lett, 2006, 97, 
H.U. Baranger, Phys. 



Liu, N. Jiang, X. Lin, W. Guo, S. X. Du, W. A. Hofer, X. 

C. Xie, and H.-J. Gao, Phys. Rev. Lett, 2007, 99, 106402 
[43] J. Lehmann, A. Gaita-Ario, E. Coronado and D. Loss, 

Nature Nanotech., 2007 2 312 - 317; J. Lehmann, A. 

Gaita-Ario, E. Coronado and D. Loss, J. Mat. Chem., 

2008 (in this same special issue, insert volume and pages). 
[44] N. Ishikawa, M. Sugita, T. Ishikawa, S. Koshihara and 

Y. Kaizu J. Am. Chem. Soc, 2003, 125, 8694-8695 
[45] M. A. AlDamen, J. M. Clemente-Juan, E. Coronado 

and A. Gaita-Arino, J. Am. Chem. Soc, 2008, DOL 

10.1021/ja801659m 
[46] N. Aliaga-Alcalde, R. S. Edwards, S. O. Hill, W. Werns- 

dorfer, K. Folting and G. Christou, J. Am. Chem. Soc, 

2004, 126, 12503-12516 
[47] R. Bagai, W. Wernsdorfer, K. A. Abboud and G. Chris- 
tou, J, Am. Chem. Soc, 2007, 129, 12918-12919 
[48] C. M. Ramsey, E. del Barco, S. HiU, S. J. Shah, C. C. 

Beedle and D. Hondrickson, Nature Physics 2008, 4 277- 

281 
[49] F. K. Larsen, E. J. L. Mclnnes, H. El Mkami, J. Over- 

gaard, S. Piligkos, G. Rajaraman, E. Rentschler, A. A. 

Smith, G. M. Smith, V. Boote, N. Jennings, G. A. Timco 

and R. E. P. Winpenny, Angew. Chem. Int. Ed. Eng., 

2003, 115 105-109 
[50] M. Affronte, F. Troiani, A. Ghirri, S. Carretta, P. Santini, 

V. Corradini, R. Schuecker, C. Muryn, G. Timco and R. 

E. Winpenny, Angew. Chem. Int. Ed. Eng., 2003, 115 

105-109 
[51] A. T. Costa, Jr., D. F. Kirwan and M. S. Ferreira, Phys. 

Rev. B, 2005, 085402 
[52] M. A. AlDamen, J. M. Clemente-Juan, E. Coronado and 

A. Gaita-Arifio, to be published 
[53] R. Raussendorf, Phys. Rev. A 2005, 052301 
[54] C. S. Lent, B. Isaksen and M. Lieberman, J. Am. Chem. 

Soc, 2003, 1056-1063 
[55] O. Waldmann, H. U. Gdel, T. L. Kelly and L. K. Thomp- 
son, Inorg. Chem., 2006, 45, 3295 
[56] A. Gaita-Arino, R Raussendorf, P.C.E. Stamp, to be 

published 
[57] M. Marcos, M. Martm-Rapiin, A. Omenat and J. L. Ser- 
rano, Chem. Soc Rev., 2007, 36, 1889-1901 
[58] Special issue of Chem. Rev., edited by C. L. Hill, 1998, 

98, 1-390 
[59] A. Miiller, P. Kogerlera and A. W. M. Dressb, Coord. 

Chem. Rev., 2001, 222, 193-218 
[60] E. Coronado, J. R. Galan-Mascaros, C. J. Gomez-Garcia, 

V. Laukhin, Nature, 2000, 408 447-449 
[61] D. Fattakhova Rohlfing and A. Kuhn, Carbon, 2006, 44, 

1942-1948 
[62] S. Caillicux, D. de Caro, L. Valade, M. Basso-Bert, 

C. Faulmann, L Malfant, H. Casellas, L. Ouahab, J. 

Fraxedas, A. Zwick, J. Mater. Chem. 2006, 13, 2931- 

2936 
[63] See for example: Y. Wang, X. Wang, C. Hu and C. Shi 

J. Mater. Chem., 2002, 12, 703 - 707 



